Theory of quantum transport in Josephson junctions with a ferromagnetic insulator 
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We investigate the Josephson transport through ferromagnetic insulators (FIs) by taking into 
account the band structure of FIs explicitly. Using the recursive Green's function method, we found 
the formation of a 7r-j unction in such systems. Moreover the atomic-scale 0-7T oscillation is induced 
by increasing the thickness of FI and its oscillation period is universal, i.e., just single atomic layer. 
Based on these results, we show that stable 7r-state can be realized in junctions based on high- 
ly superconductors with I^BaCuOs barrier. Such Fl-based Josephson junctions may become an 
element in the architecture of future quantum computers. 
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I. INTRODUCTION 

There is an increasing interest in the novel proper- 
ties of interfaces and junctions of superconductors and 
ferromagnetic materials One of the most interest- 
ing effects is the formation of a Josephson 7r-junction 
in superconductor / ferromagnetic-metal / superconductor 
(S/FM/S) heterostructures^ In the ground-state phase 
difference between two coupled superconductors is ir in- 
stead of as in the ordinary 0-j unctions. In terms of the 
Josephson relationship 

I j = I c sin 0, (1) 

where </> is the phase difference between the two supercon- 
ductor layers, a transition from the to 7r states implies 
a change in sign of Ic from positive to negative. Such a 
negative Ic was originally found in the Josephson effect 
with a spin- flip process^ - — In S/FM/S junctions, such a 
sign change of Ic is a consequence of a phase change in 
the pairing wave-function induced in the FM layer due 
to the proximity effect. The existence of the 7r-junction 
in S/FM/S systems has been confirmed in experiment by 
Ryanzanov et alZ and Kontos et al£ 

Recently, a quiet qubit consisting of a superconducting 
loop with a S/FM/S 7r-junction has been proposed^ - — 
In the quiet qubit, a quantum two-level system (qubit) 
is spontaneously generated and therefore it is expected 
to be robust to the decoherence by the fluctuation of 
the external magnetic field. From the viewpoint of the 
quantum dissipation, however, the structure of S/FM/S 
junctions is inherently identical with S/N/S junctions (N 
is a normal nonmagnetic metal). Thus a gapless quasi- 
particle excitation in the FM layer is inevitable. This 
feature gives a strong dissipative effect^— and the co- 
herence time of S/FM/S quiet qubits is bound to be very 
short. Therefore Josephson 7r junctions with a nonmetal- 
lic interlayers are highly desired for qubit application. 

On the other hand, a possibility of the 7r-junction for- 
mation in Josephson junctions through ferromagnetic in- 
sulators (FIs) have been theoretically predicted^ and in- 



tensively analyzed by use of the quasiclassical Green's 
function techniques J&il Recently, by extending these re- 
sults, we have proposed superconducting phased and flux 
qubits^— based on S/FI/S 7r-junctions. Moreover we 
have also showed that the effect of the dissipation due 
to a quasi-particle excitation on macroscopic quantum 
tunneling is negligibly small^ These results clearly in- 
dicate the advantage of the FI based 7r-junction for qubit 
applications with longer coherence time. 

However, up to now, a simple 5-function potential^ 
has been used in order to model the FI barrier. In this 
phenomenological model, the up (down) spin electrons 
tunnel through a positive (negative) delta-function bar- 
rier. Therefore, strictly speaking, this model describes 
not ferromagnetic insulators but half metals with in- 
finitesimal thickness. Moreover the possibility of the ir- 
junction formation in the finite barrier thickness case is 
also an unresolved problem. In order to resolve above 
issues, we formulate a numerical calculation method for 
the Josephson current through FIs by taking into account 
the band structure and the finite thickness of FIs ex- 
plicitly. In this paper we present our recent numerical 
results 2 ^— on the formation of the 7r-coupling for the 
Josephson junction through a FIs, e.g., L^BaCuOs and 
K2CUF4 and show that the mechanism of the 7r-junction 
in such systems is in striking contrast to the conventional 
S/FM/S junctions. 



II. MAGNETIC AND ELECTRONIC 
PROPERTIES OF FERROMAGNETIC 
INSULATORS 

In this section, we briefly describe the magnetic prop- 
erties and the electronic density of states (DOS) of FIs. 
The typical DOS of FI for each spin direction is shown 
schematically in Fig. 1. One of the representative ma- 
terial of FI is half-filled La 2 BaCu0 5 (LBCO)^- 2 ^ The 
crystal structure of LBCO has tetragonal symmetry with 
space group PA/mbm. In 1990, Mizuno et al, found that 
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FIG. 1: The density of states for each spin direction for a fer- 
romagnetic insulator, e.g., LBCO. V ex is the exchange spiriting 
and 8t is the width of the band. 
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FIG. 2: A schematic figure of a Josephson junction through 
the ferromagnetic insulators on the two-dimensional tight- 
binding lattice. 



points to a lattice site, where x and y are unit vectors in 
the x and y directions, respectively. In the y direction, 
we apply the periodic boundary condition for the number 
of lattice sites being W. 

Electronic states in a superconductor are described by 
the mean-field Hamiltonian 
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(3) 



LBCO undergoes a ferromagnetic transition at 5.2 Ki2£ 
The exchange splitting V ex is estimated to be 0.34 eV by 
a first-principle band calculation using the spin-polarized 
local density approximation^ Since the exchange split- 
ting is large and the bands are originally half-filled, the 
system becomes FI. 

An another example of FPFI is K2CUF4 compounds 
in which the two-dimensional Heisenberg ferromagnet is 
realized ! 28 ' 29 The ferromagnetic behavior of this materi- 
als has been experimentally confirmed by the magnetic 
susceptibilitj*22 and neutron diffraction measurements <2I 
Moreover a result of the first-principle band calculation^ 
indicated that K2CUF4 compounds with J aim- Teller dis- 
tortion have the electronic structure similar to Fig. 1. 
In the followings, we calculate the Josephson current 
through such FIs numerically. 

III. NUMERICAL METHOD 

In this section, we develop a numerical calculation 
method for the Josephson current of S/FI/S junctions 
based on the recursive Green's function technique.— ~— 
Let us consider a two-dimensional tight-binding model 
for the S/FI/S junction as shown in Fig. 2. The vector 

r = jx + my (2) 



Here 

K,r' = [-t s 8\ r _ r ,\^ + (-(l s + At s )5 r y\ <7 , (4) 

with 

c r = (c r f,c r ^) , (5) 

where cj. a (c r a ) is the creation (annihilation) operator 
of an electron at r with spin a = ( ~\ or i ), c means 
the transpose of c, and &q is 2 x 2 unit matrix. The 
chemical potential (l s is set to be 2t s for superconductors. 
In superconductors, the hopping integral t s is considered 
among nearest neighbor sites and we choose 

A = iAa 2 , (6) 

where A is the amplitude of the pair potential in the 
s-wave or d-wave symmetry channel, and 02 is a Pauli 
matrix. 

We consider FIs as a barrier of the Josephson junction. 
The Hamiltonian of the FI barrier is given by a single- 
band tight-binding model as 

r,r f .a r 

+ ^(4t- M + y ex )4 4 cv 4 , (7) 
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where V cyi is the exchange splitting (see Fig. 1). If V cx > 
8t (V ex < 8t), this Hamiltonian describes FI (FM). The 
chemical potential fi is set to be 



V = 



ex 

2~ 



4i. 



(8) 



The Hamiltonian is diagonalized by the Bogoliubov 
transformation and the Bogoliubov-de Gennes equation 
is numerically solved by the recursive Green function 
mcthod^r— We calculate the Matsubara Green func- 
tion in a FI region, 



9w n {r,r') f Un (r,r') 
-fL(r,r') -g*{ry) 



where 



Lo n = (2n + l)nT 



(9) 



(10) 



is the Matsubara frequency, n is an integer number, and 
T is a temperature. The Josephson current is given by 



Ij{4>) = -ietTj2 E Tr ~ Gu n (r,r')] , 

UJ n 171 — 1 

(11) 

with r' = r+x. The Matsubara Green function in Eq. © 
is a 4 x 4 matrix representing Nambu and spin spaces. 
Throughout this paper we fix T = 0.01T C , where T c is 
the superconductor transition temperature. 



IV. JOSEPHSON CURRENT FOR LOW-T c 
SUPERCONDUCTORS 

In this section we show numerical results of the Joseph- 
son current for low-T c superconductor/FI/low-T c super- 
conductor junctions and discuss the physical origin of 
the 7r-junction formation in such systems^ - — In the 
calculation, we assume t = t s and set W = 25, and 
A = A s = O.Olt. The phase diagram depending on the 
strength of V cx (0 < V ex /t < 8 for FM and V ox /t > 8 for 
FI) and Lp is shown in Fig- [3j The black (white) regime 
corresponds to the 7r- (O-)junction, i.e., 
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(12) 



In the case of FI, the 7r-junction can be formed. More- 
over, the 0-7r transition is induced by increasing the thick- 
ness of the FI barrier Lp and the period of the tran- 
sition is universal and just single atomic layer^ We 
also found that the atomic-scale 0-7r transition is also 
thermally stabled On the her hand, in the case of FM, 
the oscillation period strongly depends on V ex and the 
temperature!^ 

A physical origin of the appearance of the 7r-junction 
and the atomic scale 0-7T transition can be explained as 
follows^ In the high barrier limit (V cx ^> t), Josephson 
critical current is perturbatively given b y 20 ' 21 



I c oc T ; *T T . 



(13) 
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FIG. 3: The phase diagram depending on the strength of 14 x 
and L F for FM (0 < V ex /t < 8) and FI (F ox /£ > 8). The 
black and white regime correspond to the ir- and 0-junction, 
respectively. 



Here TWn is a transmission coefficient of the FI barrier 
for up (down) spin electrons. In the case of the single-cite 
FI (i.e., Lp = 1), the transmission coefficients are ana- 
lytically given by use of the transfer matrix method^ - — 
as 



T 



t 



oil 



Ty = 



-cti 



(14) 
(15) 



where ct\ is a spin-independent complex number. There- 
fore the sigh of the critical current 
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(16) 



becomes negative, so the 7r-junction is formed in the case 
of single-cite FI barrier. 

On the other hand, the transmission coefficients for an 
arbitrary value of Lp > 1 can be expressed by 



Ta 



OtLi 



t 



(17) 



(18) 



where oll f is a complex number. So the sign of the crit- 
ical current 



2L f 



(19) 



becomes negative for the odd number of Lp and positive 
for the even number of Lp . Therefore we can realize the 
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atomic-scale 0-7r transition with increasing the thickness 
of the FI barrier Lp as demonstrated in Fig. [3] 



boundary condition for the a and b direction and taken 
into account the <i-wave order-parameter symmetry in 
HTSC, i.e., 



V. JOSEPHSON CURRENT FOR HIGH-T C 
SUPERCONDUCTORS 

We would like to show an experimental set-up for ob- 
serving the 7r-junction using LBCO in Fig. 2J From the 
perspectives of the Fl/superconductor interface match- 
ing and the high-temperature device-operation, the us- 
age of high-T c cuprate superconductors (HTSC), e.g., 
YBa2Cu307-<5 and La2- s Sr a: Cu04(LSCO) is desirable. 
Recent development of the pulsed laser deposition tech- 
nique enable us to layer- by- layer epitaxial-growth of such 
oxide materials . 40 i 41 Therefore, the experimental obser- 
vation of the 0-7r transition by increasing the layer num- 
ber of LBCO could be possible. 

In order to show the possibility of 7r-coupling in such 
realistic HTSC junctions, we have numerically calculated 
the c-axis Joscphson critical current Ic based on a three- 
dimensional tight binding model with L a and Lb being 
the numbers of lattice sites in a and b directions [Fig. [4] 
(a)]i 23 ' 42 In the calculation we have used a hard wall 
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[cos k x a — cos k y a). 



(20) 



The tight binding parameters t and g have been deter- 
mined by fitting to the first-principle band structure cal- 
culations 2 ^. Figure 5 shows the FI thickness Lp depen- 
dence of I c at T = 0.01T C for a LSCO/LBCO/LSCO 
junction with V ex /t = 28, A^/t = 0.6, and L a = Lf, = 
100. As expected, the atomic scale 0-tt transitions can 
be realized in such oxide-based c-axis stack junctions. 

The formation of the 7r-junction can be experimentally 
detected by using a HTSC ring [see Fig. 4 (b)]. The 
phase quantization condition for the HTSC ring is given 
by 



2tt- 



27m, 



(21) 



where </>i and <p2 are the phase difference across the junc- 
tion 1 and 2, $ is the flux penetrating though the ring, 
$o is the flux quantum, and n is an integer. The current 
passed through the ring divides between the junction 1 
and 2, i.e., 



I = Ici sin 4>i + Ic2 sin </> 2 . 



(22) 



Applied external magnetic flux $ oxt depletes phases <pi 
and <p2 causing interference between currents through the 
junctions 1 and 2. For a symmetric ring with Iq\ ~ 
Ic2 = Ic an d negligible geometric inductance (L = 0), 
the total critical current as a function of $ e xt is given by 



iSP = 17? = 2Ir 
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(23) 



for the case that Lp of the both junctions are same. If 
Lp of the junction 1(2) is even and Lp of the junction 
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FIG. 4: Schematic picture of (a) c-axis stack high- 
Tc superconductor/LBCO/high-Tc superconductor Josephson 
junction and (b) high-T c ring which can be used in experimen- 
tal observations of the 7r-junction. 



FIG. 5: The Josephson critical current Ic as a function 
of the FI thickness Lf at T — 0.01T C for a c-axis stack 
LSCO/LBCO/LSCO junction with V cx /t = 28, A d /t = 0.6, 
and L a = Lb = 100. The red (blue) circles indicate the 7r(0)- 
j unction. 
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2(1) is odd, we get 



7-07T 

1 C 
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(24) 



Therefore the critical current of a 0-ir (0-0) ring has a 
minimum (maximum) in zero applied magnetic fields 
Experimentally, the half-periodic shifts in the interfer- 
ence patterns of the HTSC ring can be used as a strong 
evidence of the 7r-junction. Such a half flux quantum 
shifts have been observed in a s-wave ring made with a 
S/FM/SM and a S/quantum dot/S junction^ 

It is important to note that in the case of c-axis 
stack HTSC Joscphson junctions ; 46 ' 47 no zero-energy 
Andrccv bound-states 4 ^ which give a strong Ohmic 
dissipation 4 ^— are formed. Moreover, the harmful in- 
fluence of nodal-quasiparticles due to the d-wave order- 
parameter symmetry on the macroscopic quantum dy- 
namics in such c-axis junctions is found to be week both 
theoretically^— and experimentally^ 7 - - — Therefore 
HTSC/LBCO/HTSC 7r-junctions would be a good can- 
didate for quiet qubits. 



VI. SUMMARY 

To summarize, we have studied the Josephson effect in 
S/FI/S junction by use of the recursive Green's function 
method. We found that the 7r-junction and the atomic 
scale 0-7r transition is realized in such systems. By use 



of the transfer matrix calculation, the origin of the 7r- 
j unction formation can be attributed to the n phase dif- 
ference of the spin-dependent transmission coefficient for 
the FI barrier. Such FI based 7r-junctions may become 
an element in the architecture of quiet qubits. 

We would like to point out that the 7r-junction can be 
also realized in the Josephson junction through an an- 
other type of FI, i.e., a spin- filter material, in the case of 
the strong hybridization between localized and conduc- 
tion electrons . 61 ' 62 It should be also note that FI mate- 
rials treated in this paper can be categorized in strongly 
correlated systems. Moreover, in actual junctions, the 
influence of the interface roughness could be important. 
Therefore investigation of the atomic-scale 0-7r transition 
in the presence of the many-body and disorder effect will 
be also the subject of future studies. 
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